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ABSTRACT 



We present the longest, by a factor of two, near-infrared lightcurve from 



<N ; Sgr A* - the supermassive black hole in the Galactic center. Achieved by com- 

^ i bining Keck and VLT data from one common night, which fortuitously had si- 

OO ! multaneous Chandra and SMA data, this lightcurve is used to address two out- 

^ ' standing problems. First, a putative quasi-periodicity of ~20min reported by 

groups using ESO's VLT is not confirmed by Keck observations. Second, while 
^ I the infrared and mm-regimes are thought to be related based on reported time 

\ lags between lightcurves from the two wavelength domains, the reported time lag 

of 20min inferred using the Keck data of this common VLT/Keck night only is 
at odds with the lag of ~ 100 min reported earlier. With our long lightcurve, we 
find that (i) the simultaneous 1.3 millimeter observations are in fact consistent 
with a ~ 100 min time lag, (ii) the different methods of NIR photometry used 
by the VLT and Keck groups lead to consistent results, (iii) the Lomb-Scargle 
periodogram of the whole NIR lightcurve is featureless and follows a power-law 
with slope -1.6, and (iv) scanning the lightcurve with a sliding window to look 
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for a transient QPO phenomenon reveals for a certain part of the hghtcurve a 
25 min peak in the periodogram. Using Monte Carlo simulations and taking the 
number of trials into account, we find it to be insignificant. 

Subject headings: black hole physics, Galaxy: center 



Introduction 



The near-infrared (NIR) regime at high angular resolution has proven to be of great value 
in Galactic center research. The proper motion of stars detected in this waveband demon- 
strates the existence of a supermassive black hole (BH) at the center of our Galaxy: Sagittar- 
ius A* (Sgr A*: see, e^ Eckart fc Genzell EqQgI : iGhez et ahlEoOsl . boool . I2OO5I : iGenzel et al. 

20021). In 2003 . NIR emission associated with Sgr A* was detected 



2000 



Schode 



et a 



( IGenzel et al.ll2003l : iGhez et al.ll2004l ). which is important since it is the most underluminous 
BH accretion system observed thus far (with a bolometric luminosity nine orders of magni- 
tude lower than its Eddington luminosity). The NIR emission is highly variable: intensity 
changes by factors < 10, lasting between about 10 an d 100 min, occur a t least 4 times a 
day, and are highly polarized (e.g. Eckart et al. 2006a| jbl: iMeyer etaL 2006al lbl: IDo et al. 2008 



Trippe et al.ll2007l ). Simultaneous NIR and X-ray observations have revea led that each X-ray 
flare is accompanied by a NIR flar e with zero time lag, bu t not vice versa feckart et al. 2004, 
2006al . I2OO8I : iBelanger et al.l l2005l : iHornstein et al.l 120071 ) . Recent campaigns that included 
mm observations, reported a characteristic time lag between the NIR/X-ray and longer 
wavelengths that has been interpreted in terms of an expanding synchrotron plasmon. The 
limited overlap between data sets, however, has led to debates over the nature of this time 
lag (e.g. lEckart et al.1l2006al : IVusef-Zadeh et al.lbood . I2OO7I : iMarrone et aPboOTh . 



The characteristics of the NIR emission from Sgr A* is of particular interest given 

^20 min present 



Genzel et al. 


2003; 


Eckart et al. 


2006b 



robust statistical estimators and found no signiflcant peaks in the periodograms of the Sgr A* 
lightcurves observed with the Ke ck II telescope. Moreover , quasi-periodicities in the X-ray 



regime claimed to be present by lAschenbach et al.l (120041 ). are not statistically signiflcant 



(IBelanger et al.ll2008[ ). 



In this Letter, we address the questions of a periodic component in the NIR flux, and the 
relation of the NIR to the mm-regime, by combining for the flrst time contiguous Keck and 
VLT data. During that 10 hour session, there was also simultaneous coverage by Chandra and 
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SMA. These observations were published by iMarrone et aLl (120071 ). but interpreted taking 
only the Keck data into account. 



The data 



During the night of 2005 July 30-31, the VLT observed Sgr A* from 23:05 UT to 06:53 UT 
using the Natural Guide Star Adaptive Optics system and NIR camera NACO on 
Since these observations have not been published before, the observational details are given 
here. The detector integration time was 15 s, and four images were co-added before the data 
were recorded: the effective resolution is ~1 image/80 s. Dithering was used to minimize 
the effects of dead pixels. The atmospheric seeing conditions ranged between 1" — 2.5" 
(as determined by the differential image motion monitor) during the first two hours, and 
0.5" — 1.25" afterwards (the average Strehl ratio of the VLT data is 17%). 

At 07:00 UT, the Keck II telescope started its monitoring campaign using the NIRC2 
camera in combination with Laser Guide Star A O. The Keck observat ions were made by 
cycling through the H-K-L wavelength filters (see iHornstein et al.l 120071 . for details). Since 
the VLT observations were carried out in K-band, we consider only the K-band subset of 
the Keck data. 

Both data sets were reduced in the same standard way, i.e. sky subtracted, flat-fielded, 
and corrected for bad pixels. Images with a Strehl ratio less than 10% were removed (13 
out of 266 VLT images). For every i ndividual imag e , the point spread function (PSF) was 
extracted with the code StarFinder by lDiolaiti et al.l (120001 ). Each exposure was deconvolved 
with a Lucy-Richard deconvolution and restored with a Gaussian beam. The flux of Sgr 
A* and other compact sources in the field were obtained via aperture photometry on the 
diffraction limited images with a circular aperture of radius 0".03. The background flux 
density was determined as the mean flux measured with apertures of the same size at five 
different positions in a field located 1" northwest of Sgr A* that shows no individual stars. 
Photometric calibration was done relative to stars in the field with known flux. For the 
extinction correction we assumed Ak = 2.8 mag (lEisenhauer et al.l 120051 ). Estimates of 
uncertainties were obtained from the standard deviation of fluxes of nearby constant sources. 

It is noteworthy that this is the first time that Galactic Center VLT and Keck data have 
been reduced homogeneously. While the VLT groups mainly use deconvolution and aperture 
photometry, the Keck group uses PSF fitting without deconvolution. The Keck part of the 
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lightcurve presented in the next section is similar to that reported by lHornstein et al.l (120071 ). 
where the PSF fitting technique was apphed. After rescahng the lightcurve by a multiplica- 
tive factor to account for a dif f erent de-reddening factor and slightly different calibration 
values used by iHornstein et al.l (120071 ). the average difference per data point between both 
methods is only 0.031 mJy (de- reddened), well within our Icr-error bar of 0.18 mJy. Both 
data reduction methods are consistent with each other, and no large systematic errors are 
introduced by choosing either one. 



3. Results 



3.1. The NIR properties of Sgr 



Figure [T] (left panel) shows the 2005 July 30-31 de-reddened lightcurve. The first 450 
min are the VLT data, and the following 130 min are the Keck data. Note that the error 
bars are smaller for the Keck data due to better seeing conditions at Mauna Kea that 
night. The corresponding Lomb-Scargle periodogram for the ent ire lightcurve is presented 
in the right panel of Fig. [U It is consistent with the finding of iDo et al.l (120081 ) that the 
NIR emission of Sgr A* is described by a single stochastic process (other than additive 
measurement noise) that has a power-law spectrum, V cx /~", very similar to the X-ray 
emission of AGN. The Lomb-Scargle periodogram in Fig. [1] is a single realization of this 
process and therefore fluctuates around the spectrum. We determine the probability density 
function (PDF) of these fluctuations around the spectrum at a given frequency empirically 
with Monte Carlo simulations. With the PDF at hand, we can assess the likelihood that 
peaks in the periodogram are not due to a fluctuation but rather have a physical cause 
intrinsic to the source. 



We used two di f ferent Monte Carlo based analyses recently developed by lDo et al.l (120081 ) 
and iBelanger et al.l (120081 ) to look for significant peaks in the periodogram. These analyses 
are carried out by first determining the power-law index of the sp ectrum of the stochastic 
process, and then generating lightcurves (realizations of this process; iTimmer fc Koniglll995l ) 
with a sampling function matching those of the data set. Finally, the significance of each 
observed periodogram peak is derived from the large simulated reference data set. 

An accurate determination of the power-law i ndex a of the underlying stochastic process 
is important in this approach. iBelanger et al.l (120081 ) do this by first performing several 
estimates of a by fitting a power-law to the periodogram made from the lightcurve binned 
with successively larger bin times, and then carrying out simulati ons with t he sam e count 
rate and sampling to find the matching curve of a versus bin time. iDo et al.l (120081 ) also use 
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MC simulations, but determine the power-law index using the structure function, which also 
takes the sampling into account. Both methods lead to a spectral index of a = 1.6 ± 0.05 
(formal fitting error). Note that the length of the lightcurve allows to sample frequencies 
< 10~^ Hz (0.006 min~^) for the first time, showing that the power-law extends to this regime. 
Unfortunately, we cannot fully exclude the possibility of spectral leakage from a process with 
a spectral index steeper than -1.6. However, our main conclusions in this paper also hold 
true for steeper power-law indices. We plan to use more sophisticated spectral estimators 
elsewhere. 

The periodogram in the right panel of Fig. [1] clearly shows no outlying peaks above 
the underlying spectrum of the power-law process. However, the featurelessness of the peri- 
odogram of the whole lightcurve does not rule out the presence of a periodic component in 
parts of it, as the mechanism giving rise to such a component may be transient and short 
lived. We therefore did a search for periodicities over a restricted range of frequencies by 
scanning the lightcurve using a sliding window method where window here means a sub-span 
of the time series. This consists of constructing a periodogram for the data subset corre- 
sponding to each window, and assigning a significance to each point based on the probability 
density functions derived from the simulations of the whole lightcurve described above. 

Fig. [2] shows the result of such a scan using a 60 min window with 5 min steps. The prob- 
ability that the most significant peak in each periodogram is due to a statistical fluctuation 
around the power-law spectrum is plotted against the start time of the sliding window. The 
most significant peak overall (the one with the lowest probability to be due to a fiuctuation) 
is found in the window starting at minute 385, and we thus looked at this window subset in 
more detail. Indeed, the fiu x between 385 and 4 4 5 min looks v e ry sirn ilar to t he 'sub-flare' 
phenomenology reported by iGenzel et al.l (120031 ). lEckart et al. (EoSQ), and iMeyer et al.l 



( l2006al ). These sub-flares are flux peaks superimposed on broader, longer lastin g flux excur 



sions a nd are t hought to be the rn anifestation of the claimed quasi-periodicity. iMeyer et al. 



fl2006al lb[) and lEckart et al.l (120081 ) showed that they can be interpreted in terms of a rel- 



ativistically orbiting spot whose emission adds to the emission of the accretion flow. The 
NIR flux can therefore be described as F(t) = A{t) + M(t) + S{t) with A{t) a stochastic 
process with a power-law spectrum as above, M{t) uncorrelated measurement noise, and 
S{t) the deterministic flux of an (evolving) orbiting spot leading to a (quasi-)periodicity. 
Our MC simulations include A{t) and M(t) so that a possible quasi-periodic component can 
be identifled as signiflcant in the periodogram. 

Fig. [2] shows that there is a peak in the periodogram corresponding to the 385 — 445 min 
subset which has a false alarm probability of only 2 • 10~^ (this corresponds to A.2a in 
Gaussian equivalent terms, but note that the PDF is not Gaussian). This peak occurs at a 
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frequency of 0.04 min~^ (25min). However, we have to ask the question how hkely it is that 
a periodogram peak with such a low probabihty does not only occur at this window, but in 
any window. After all, there is a priori nothing special about this certain sub-span of the 
lightcurve. We therefore count how often a peak with probability < 2 • 10~^ occurs in our 
simulations in any window: for a fixed window length of 60min, we find 3120 occurrences 
for 30,000 simulations. This implies that the overall false alarm probability of the peak in 
the 385 — 445 min subset is 0.104 (corresponding to 1.6a). Furthermore, we have to take into 
account the trials with windows of different lengths in scanning the data: this immediately 
yields a final significance of < la. 

Hence, we conclude that the whole lightcurve is consistent with a pure power-law process 
and no periodic component is needed. It is important to point out, however, that if we had 
observed (or analyzed) only the part of the data between 385 — 445 min, our result would 
have been interpreted as a A a detection of a 25 min QPO. Thi s point s to an explanation 
of the c ontradictory results of 
( l2006al ) , on t he one hand, and 



Genzel et all J2003h . Ifickart et all J2006al lb[) and iMever et al. 



DoetaL 



(120081 ). on the other. Our results here seem to favor 



the finding of iDo et al.l (l2008l ) where no periodicity was detected. 



It is interesting to note that the 25 min periodogram peak is most significant for a 
60 min sliding window, implying that only the first two of the four 'sub-flares' between 380 
and 460 min are sampled. The reason for this is probably that the period - if present - 
is evolving over the four cycles resulting in a periodogram peak which is too wide to be 
identified as a significant periodic component. This, however, does not exclude the presence 
of an evolving, inwards spiraling bright sp ot or more complica ted hydrodynamic instabilities 
in the accretion flow around Sgr A* (e.g. iFalanga et al.ll2007l ). 



3.2. The connection to the miUimeter regime 



The Chandra X-ray Observat ory and the SMA obse rved Sgr A* simultaneously to this 
new NIR lightcurve as reported by lHornstein et al.l (120071 ) and lMarrone et al.l (120071 ). Chan- 
dra started its observations at 20:00 UT (30 July) and stopped at 08:30 UT (31 July), so that 
there are simultaneous X-ray data for the entire NIR lightcurve seen in F ig. [H No X-ray 
f lare o ccurred during the overlap time (23:05 UT - 09:00 UT), see Fig. 6 in iHornstein et al. 
( 120071 ). despite the activity in the NIR with flares which are somewhat stronger than aver- 
age. This may be due to the higher X-ray background caused by the steady Bondi-Hoyle 
accretion flow within 1" around Sgr A* or/and a large surface area of the flaring r egion that 
leads to a low synchrotron self-Compton luminosity (see also lMarrone et al.l 120071 ). 
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In contrast to the X-ray lightcurve, the SMA 1.3 mm observations, which started at 
05:28 UT, show some variabihty (see Fig. [3]). In particular, they show one fla re at 08:20 UT, 



rough ly 20min after the NIR flare observed with Keck. In their analysis, iMarrone et al. 



( 120071 ) took only the Keck NIR data into account. They interpreted this 20 r ain time delay 
i n term s of an expansion of energetic plasma similar to the model proposed by Ivan der Laan 
( 119661 ). In such a model, the expanding plasma region is optically thick before, and optically 
thin after the peak flux. As this transition is frequency dependent, smaller and later flare 
peaks are expected at longer wavelengths. The 20 min la g, however, diff e rs fro m other 
reported time lags. In a second data set from 2006 July, iMarrone et al.l (120071 ) infer a 
97 ± 17 r ain time lag between an X -ray flare (which is sup posed to be synchro i ious w ith the 
NIR, see lEckart et al.ll2006al . l2008l ) and a millimeter flare. lYusef-Zadeh et al.l (120071 ) infer a 
lag of 110 ± 17 min for the same data. 

With the full NIR information at hand, the assumption of a correlation between the 
08:00 UT NIR flare and the 08:20 UT millimeter flare becomes uncertain. Fig. [3] shows that 
there is another equally bright NIR flare preceding the 8:00 UT flare, with no counterpart 
at a 20 min lag in the 1.3 mm lightcurve. Therefore, it is equally likely that the mm flare is 
correlated with t he preceding flare, as it is with the 8:00 UT flare: substantially altering the 
interpretation of IMarrone et al.l (120071 ). 



In fact, the huge difference in the reported time lags in the 2005 mm-IR and the 2006 
mm-X-ray data can be resolved with the NIR data presented here: if the ~8:20UT mm 
flare (Fig. [3]) is correlated with the wide NIR flare at ~6:00UT, then the inferred time lag 
is ~140min. This la g is much closer to 97min/110 min, than to the 20 min assumed by 
Marrone et al.l (120071 ). Note that the 8:00 UT NIR flare might then be correlated with the 
SMA flare seen at 10:15 UT. 



4. Conclusions 

In this Letter we report on a 600 minute NIR lightcurve of Sgr A* that combines Keck 
II and VLT data. We showed that the Lomb-Scargle periodogram of the overall lightcurve is 
featureless, and is statistically consistent with a single stochastic process that has a power- 
law spectrum of index —1.6. A certain subset of the data has a prominent peak in its 
corresponding periodogram which, however, is not signiflcant when analyzed in the context 
of the whole lightcurve. This points to the following dilemma in Sgr A* research: if a periodic 
component exists, it is clearly a weak and transient phenomenon that persists over very few 
cycles, and probably has an evolving period. As we have shown, certain parts of longer pure 
red noise lightcurves can easily mimic such a behavior. To be able to distinguish between 
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both scenarios, many lightcurves are needed and the range of frequencies under consideration 
must be narrowed down, e.g. to 0.04 — 0.07 min~^ if a noticeable peak continues to occur in 
this range only. 

We furthermore showed the difficulty with establishing the simple expanding plasmon 
model frequently proposed to relate NIR and mm-flares. Sgr A* is such an active source 
in the NIR (and maybe not every NIR flare has a millimeter counterpart) that very long 
simultaneous observations are needed for a meaningful cross-correlation analysis. 

We are very grateful to Dan P. Marronc for providing us with the SMA data. Some of 
the data presented here were obtained from Mauna Kea observatories. We are grateful to 
the Hawai'ian people for permitting us to study the universe from this sacred summit. This 
work was supported by NSF grant AST-0406816. 

Facilities: VLT:Yepun (NACO), Keck:II (NIRC2), SMA 
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Fig. 1. — Left panel: De-reddened lightcurve from Sgr A* (lower curve) on the night of 2005 
July 30-31 (UT) and from a comparison star (SO-8 shifted 5mJy upwards; upper curve). 
The circles are the VLT data points (0-450 min) and the triangles are the Keck data points 
(450-580 min). The small gaps in the hghtcurve are due to either periods of very bad seeing 
leading to useless data, sky observations, or laser collision at Keck. Right panel: Lomb- 
Scargle periodogram of the lightcurve. No smoothing was applied. 
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Fig. 2. — Sliding window scan of the lightcurve seen in Fig. [TJ Each time window has a length 
of 60 min and has been shifted by 5 min for the scan. The abscissa shows the start time of 
each sliding window. The ordinate shows the probability that the most significant peak in the 
corresponding periodogram is due to a fluctuation around the power-law spectrum of the red 
noise stochastic process. A very low probability means that an additional (quasi-)periodic 
component is needed to explain the periodogram peak. See text for more details. 
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Fig. 3. — The SMA 1.3 mm data (upper curve) from iMarrone et al.l (120071 ) plotted over a 
common time axis with the overlapping part of the NIR lightcurve (lower curve) shown in 
Fig. [H The flux units are mJy for the NIR and Jy for the SMA data which have been shifted 
upwards by 3.5 Jy for better comparison. Please note that the first four points in the SMA 
lightcurve are somewhat unreliable (D. Marrone, priv. comm.). The scatter in these four 
points seems large compared to the error bars, which indicates that they should in fact be 
larger for these points than shown in the plot. 



